Goal:

Questions:

What does it mean to be correlated? Might this occur if I simulatate some controlled movement?

Background

Experiment

I created a simulation study to observe the co-evolutionary outcome of the newt-snake interaction with different genetic architectures (GAs) in a spatial setting. I hypothesized that we would see an interaction (co-evolutionary arms race) between newt and snake phenotype under some GA combinations when newts and snakes were evolving over geographical space. Each GA is paired with another GA creating 16 combinations.

GA1 experiment values:

Each GA combination and trial has its own msprime simulation, but the msprime file is shared between the 5 section and 7 section slim runs.

The data

## All cor, lit, and grid files exist!
## This program will now end!

Mean Phenotype Whole Simulation

Phenotype differences

Table of average Differences

##                    Group.1           x
## 1  1e-08_0.005_1e-08_0.005 -0.03323921
## 2   1e-08_0.005_1e-09_0.05 -1.30511227
## 3    1e-08_0.005_1e-10_0.5 -2.50095446
## 4      1e-08_0.005_1e-11_5 -1.19918146
## 5   1e-09_0.05_1e-08_0.005  1.33051637
## 6    1e-09_0.05_1e-09_0.05 -0.43179017
## 7     1e-09_0.05_1e-10_0.5 -1.13158335
## 8       1e-09_0.05_1e-11_5 -0.48635858
## 9    1e-10_0.5_1e-08_0.005  2.03584772
## 10    1e-10_0.5_1e-09_0.05  0.39730616
## 11     1e-10_0.5_1e-10_0.5 -0.73000969
## 12       1e-10_0.5_1e-11_5 -0.40944406
## 13     1e-11_5_1e-08_0.005  0.98550774
## 14      1e-11_5_1e-09_0.05 -0.79978450
## 15       1e-11_5_1e-10_0.5 -1.42371425
## 16         1e-11_5_1e-11_5 -0.60451822

Connection between higher phenotype and population

Phenotype differences

Phenotype & Populationsize differences

Correlation

cor

Correlation Histograms

Plot 1

Plot 2

Plot 3

Plot 4

Plot 5

Plot 6

Plot 7

Plot 8

Plot 9

Plot 10

Correlation across time

Random 1

## [1] "pattern 1e-10_0.5_1e-08_0.005_3"
## [1] "Cor between average snake pheno and local cor -0.15534739530592"
## [1] "Cor between average newt pheno and local cor 0.141016273773035"
## [1] "Cor between average dif pheno and local cor -0.253303256208326"
## [1] "Cor between newt pheno and snake 0.746100350522642"

Random 2

## [1] "pattern 1e-10_0.5_1e-09_0.05_2"
## [1] "Cor between average snake pheno and local cor 0.434830748796898"
## [1] "Cor between average newt pheno and local cor 0.303736141742003"
## [1] "Cor between average dif pheno and local cor 0.121715259013478"
## [1] "Cor between newt pheno and snake 0.897084082842825"

Random 3

## [1] "pattern 1e-08_0.005_1e-08_0.005_3"
## [1] "Cor between average snake pheno and local cor -0.423984702204067"
## [1] "Cor between average newt pheno and local cor -0.262198745181791"
## [1] "Cor between average dif pheno and local cor 0.200962706000563"
## [1] "Cor between newt pheno and snake 0.688195150229443"

What happens over time (looking at the beginning, middle, and late part of my simulations)

Pheno Beginning

Pheno Middle

Pheno End

Dif Beginning

Dif Middle

Dif End

Summary

Early-Sim Population Size Summary

Mid-Sim Population Size Summary

Late-Sim Population Size Summary

Early Difference Summary

Mid Difference Summary

Late Difference Summary

By Snake GA (Early)

By Snake GA (Mid)

By Snake GA (Late)

By Newt GA (Early)

By Newt GA (Mid)

By Newt GA (Late)

Heatmap

Population Size (Early)

Population Size (Mid)

Population Size (Late)

Phenotype (Early)

Phenotype (Mid)

Phenotype (Late)

What is up with the correlations

Early Simulation Correlation

Mean

## [1] 0.4002541

Max

## [1] 0.4181679

Min

## [1] 0.1298847

popsize

## [1] 0.5833832

Late Simulation Correlation

Mean

## [1] 0.200901

Max

## [1] 0.1563885

Min

## [1] 0.09847233

Pop size

## [1] 0.409964